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Abstract 

Let {D{s),s > 0} be a Levy subordinator, that is, a non-decreasing process with stationary and 
independent increments and suppose that D{0) — 0. We study the first-hitting time of the process D, 
namely, the process E{t) = inf{s : D{s) > t}, t > 0. The process E is, in general, non-Markovian 
with non-stationary and non-independent increments. We derive a partial differential equation for the 
Laplace transform of the n-time tail distribution function P[E{ti) > si, . . . , E{t„) > s„], and show that 
this PDE has a unique solution given natural boundary conditions. This PDE can be used to derive all 
n-time moments of the process E. 



1 Introduction 

Consider a non-decreasing Levy Process {D{s), s > 0}, starting from 0, which is continuous from the right 
with left hmits. Such a process is called a subordinator. It has stationary and independent increments and 
is characterized by its Laplace Transform 

The function (f> is called the Laplace exponent and is given by the Levy-Khintchinc formula: 

0(A) = AiA+ [ {l-e-^^)U{dx), (1) 

J(0,oo) 

where /i > is the drift and H is a measure on IR+ U {0} which satisfies /q°°(1 A x)Il{dx) < oo (see [T], [7] 
or 0). 
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We want to study the first-passage time of such a process, focusing on its finite-dimensional distributions. 
The first-passage time of a subordinator {D{s),s > 0}, is a new process {E{t),t > 0}, commonly called an 
inverse subordinator, and is defined as follows: 

E{t) = inf{s : D{s) > t}, t> 0. 

It is worth noting that in certain cases, the Levy subordinator D is itself a first passage time of another Levy 
processes. For example, the ^-stable subordinator is the first passage time of standard Brownian motion, 
and the inverse Gaussian subordinator is the first passage time of standard Brownian motion with drift ([1], 
exercise 2.2.10). For these examples, E is the "first-passage time of a first-passage time". 

Inverse subordinators appear in a variety of applications. For instance, they are used in the study of 
scaling limits of continuous-time random walks and fractional kinetics, U, [S], [T3], [13 > El- Here, the time 
variable of a Markov process (typically Brownian motion), is replaced by an inverse a-stable subordinator, 
with < a < 1. This particular time change gives rise to anomalous diffusion, or sub-diffusion, where the 
variance of the process grows at a rate which is non-linear in time. The a-stable subordinator is "fast" 
because all moments of order a or less are infinite. An "ultrafast" subordinator introduced in [16 and [13] 
is a process which has, in general, no finite moments. While these subordinators have infinite moments, 
the moments of their inverses are finite. For the a-stable subordinator, the first moment of its inverse 
subordinator, EE{t), grows as t". For the ultrafast subordinator, the inverse is "ultraslow", and can grow 
like a slowly varying function. 

Our goal in this paper is to characterize the n-point distribution function of a general inverse subordinator 
with a simple partial differential equation and to use this PDE to derive explicit expressions for joint moments. 
A PDE was derived (heuristically) by [Ij, [S] for the joint density function in the special case of a-stable 
inverse subordinators and was used to obtain joint moments of these processes. To obtain precise results in 
the general case, we found it convenient to use suitable modified cumulative distribution functions instead 
of densities. An alternative approach, using Cox processes was followed by |15j . who also obtains joint 
moments of increments of inverse subordinators. Our results are implemented in |19j . where we develop 
algorithms for computing the moments numerically. 

This paper is organized as follows: The n-point distribution of E is studied in Section [2] where we derive 
a PDE for the Laplace transform of the n-point tail distribution function P[E{ti) > si, . . . ,E{tn) > s„]. 
In Section [3] we use this PDE to calculate the n-point moments of a general inverse subordinator. The 
examples presented briefiy in Section |4] are discussed in detail in the companion paper |19j . For convenience, 
we provide an appendix which contains key path properties of inverse subordinators. 

2 PDE for the Multiple-time Distribution Function of Inverse 
Subordinators 

Unlike the Levy process D, the inverse subordinator E is often non-Markovian. Thus, we must consider all 
finite-dimensional distributions. For n > 1, define 

il(")(si,...,s„,ii,...,i„) = P[^(t,) >s„ z=l,...,n]. (2) 

To simplify notation, we will write _ff(")({si, tijf^i) = ff(")(si, . . . , s„, ti, . . . , t„). In the following, M" = 
{(si, . . . , s„), Si > 0,i = 1,...,??,}. Also, C^(M!J:) will denote the space of continuously differentiable 
functions on M" . The following theorem shows the Laplace transform, 

poo poo 

^("H{s^,Aar=i)= / ••■/ e-^^'^--^-'-H^^H{s,,U}7=i)dti...dt. 
Jo Jo 
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of satisfies a PDE. Observe that HW{0,X) = J^' P[E{t) > Oje'^^dt = 1/A since P[E{t) > 0] = 1 for 
t > 0. We let H^'^\{si = 0}) denote ij'"^ in ^ with = and the other arguments unchanged. 

Theorem 2.1 Let D be a general Levy suhordinator and let E he the inverse subordinator of D. For 
Si,S2,---,Sn > 0, the Laplace Transform of the n-point tail distribution iJ*^"^ of E defined by ^ is the 
unique solution in C^(K") to the following PDE 

^— + TT- + • • • + TT") ^ = -'^(^1 + • • • + >^n)H^\ (3) 

^dsi ds2 OSnJ 

together with the boundary conditions 



i7(i)(0,A) = i, n = l, 

H^H{s^^O})^ J-H^^K{sk,Xk}k^^), z=l,...,n, n>l. 

Proof. From Proposition |A.2[ we have that 

H^''\{s^,t,}^^,) = P[E{t,)>s,, t^l,...,n]. 

= P[D{si) <t.,, i^l,...,n], for a.e. ti, . . . ,t„ > 0. 

Let So = and fix si,S2, ■ ■ ■ , Sn > 0. Calculating the Laplace transform in ti, . . . , t„ of we have 

^ poo poo ^ 

Jo Jo 

/ 1 \n poo i-oo I Qn " \ 

Integrating by parts in ti, . . . , we obtain 

— (_l)2n 



(4) 



) < tj, i = 1, . . . , . . . dtn- 



^ Eexp(-^A,£>(s,)). 



Now rewrite the variables so,si,...,s„ is an increasing order, = Sj(o) < Sj(i) < ■•• < ^j{n)^ where 
j(0) = and j(l), . . . tjin) is a permutation of the integers 1, . . . , n. We add and subtract terms in the 
summation above and use the independence and stationarity of the increments of D to obtain 



i/(")({s„A,}r=i) = j^^Eexp (-f] (f^A,(fe)] (D(s,(,)) 



i— 1 \k—i 



1 " / / " 



(5) 
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Differentiating with respect to Si, i — 1, ... ,n, we have that satisfies the system of PDEs 

k=i+l 



"''^^(^ \ k=i 



Adding these n equations together gives 



Since j is simply a permutation of the integers {1, . . . , n}, this is equivalent to^^. 

Now we claim that with the boundary conditions (^4|, the tail probability is the unique solution to 



the PDE pi in C^{Wl). Suppose Hl{{s^, X,}f^^) is another solution in C^{Wl) to ([3| with the boundary 

conditions (Ij). Define V{{s^, X^}^^^) = - i/i)({sj, AJf^J. By linearity, V is also a solution to ^ 

with boundary conditions 

V{{s,^0}) = 0, i^l,...,n. (6) 

Fix (si,...,s„) € and WLOG, assume si < ■ ■ ■ < (if not we can simply re-index as before). 
Proceeding by the method of characteristics ([12], section 3.2), we define the function 

^.(T) = V^({s,+T,Ajr=l) T>-Sl. (7) 

Observe that (|6| can now be written as 

y(0,S2,...,s„,Ai,...,A„) = 0. (8) 
Differentiating (|7| with respect to r and using the definition of we have 



^dsi ds2 ds. 

Since V satisfies the PDE (|3|, equation ([9| implies that v satisfies the ODE 

z/'(t) = -<?!)(Ai + --- + A„)i^(t). 

Now, using r = — si as a initial condition and using (|8|, we have 

iy{-si) = V{0, S2 - si, . . . , s„ - si, Ai, . . . , A„) = 0. 

Thus, i^(r) = for r > — si is a solution to this initial value problem, and by standard uniqueness theory of 
ODEs ([9J, chapter 2), this is the unique solution. Hence 

= i.(0) = V{{s,, A,}r=i) = S>)({s„ Aar=i) - H,i{s,, AJILi), 
implying is the unique solution to with the boundary conditions ml. 
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Remark. The boundary conditions Q are natural ones, since they imply that the n-dimensional distribution 
function can be reduced to the n — I dimensional distribution function when Si = for some i. 

Remark. One might wonder why a PDE for the n-time Laplace transform of iJ*^"^ is useful when in fact 
we can write down its solution ([s]) in closed form. It is, because: 

• It simplifies the calculation of the moments of E, which we demonstrate in Section [3j 

• Understanding the dynamics of E is useful in the study of more complicated processes. For example, 
in [S], the PDE ([s]) in the case of the a-stable process (where 0(A) = A") is used to derive equations 
corresponding to the n-time distribution functions of the so-called fractional kinetic process Za(t) — 
B{E{t)), where B is Brownian motion and E is an inverse a-stable subordinator. This process appears 
as a scaling limits for various trap models, [2J. Thus, the PDE ([3| can extend this analysis done in [5] 
to a larger class of processes. 



3 Moments of Inverse Subordinators 



In this section we use Theorem 2.1 to calculate moments of a general inverse subordinator. The utility of 
the PDE ([3]) will become clear in that it will simplify many of the following computations. 

Before we calculate moments, we first argue that all moments of an inverse subordinator are finite. Notice 



that, for any a; > 0, we can bound the tail distribution of E using (31) from the proof of Proposition A. 2 
and Markov's inequality: 

P[E{t) > s]< P[D{s) <t] = P[e-^^(") > e"^*] < e^*Ee-^^(") = e^*e-"'^(^\ (10) 

which implies that EE{t)'' < 00 for any 7 > 0. 

The Laplace transform of KE{ty' with 7 > has a simple form: 



EE{ty 



poo poo 

' / e-^* / s'^-^P[E{t) > s]dsdt 
Jo Jo 

/>oo />oo 

' / s'^-i / e-^\l - P[E{t) < s])dtds 
Jo Jo 



1 I s-'-'e-'^^'Us 







7r(i + 7) 

A0(A)-r ' 



(11) 



(12) 



where we have used integration by parts to calculate e '*'*(1 — P[E{t) < s])dt. Of particular importance 



is the mean of E(t). Let U{t) — E£'(i). From (12) U has Laplace transform given by 



We will see in the following that U characterizes all finite-dimensional distributions of the process E. While 
the Laplace transform of U is easy to express in terms of 0, calculating the inverse is, in general, not always 
an easy task. In [19 , a numerical method is given for computing U for general (jj. 
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Remark. Although the subordinator D is, in general, different from a renewal process, it has some of its 
characteristics. For example, it satisfies the so called renewal theorem, which states that if the mean of the 
subordinator is finite, then ¥.E{t) = U{t) ~ WHIT)^ as t ^ co. This was proven using various methods in 
[15] . [8] and [11]. This fact can be easily seen from the Laplace transform of U given above. Indeed, if 
-D(l) has finite mean, then from the Levy-Khintchine formula, 

ED{1) = (/)'(0) = lim = + / x\l{dx) < oo, (14) 
A^o A Jq 



and thus 0(A) - ED{1)\ as A ^ 0. From ([l3|, U{\) ~ 1/{ED{1)X'^) as A ^ and the Tauberian theorem 
([7], page 10) implies the renewal theorem: 

m-^y t^OO. (15) 

Thus, for subordinators with finite mean, their mean first-passage time will exhibit a non-linear transient 
behavior for small times, followed by a linear behavior for large times. 



As an application of the differential equation given in Theorem |2.1| we obtain expressions of the Laplace 
transforms for the n-time integer moments of E in terms of Laplace transforms of lower-order moments. 

Theorem 3.1 Let D be a general Levy subordinator with Levy exponent cf) and let E be the inverse subor- 
dinator of D. For positive integers mi, . . . , to„, let 

U{ti, . . . , mi, . . . , m„) = EEit^r^ . . . E{t^r- . (16) 

Ln the special case n = I, mi = I, we will simply write U(t,l) = U{t). The n-time Laplace Transform of U 
is given in terms of strictly lower order moments by 

1 " _ 

[/(Ai, . . . , A„; mi, . . . , m„) — - — ; — — miU{Xi, . . . , A„; mi, . . . , mi_i,mi - 1, m^+i, . . . , m„). 

(p{\i H h A„) ^ 

(17) 

Proof. We first write U in terms of the tail probability H^'^\ 



oo 



C/(ti,...,t„;mi,...,m„) = Af / ... / s^-^ . . . ^ J-U)dsi, ■ • ■ , ds„, 

Jo JQ 

where M = Y[i=i "^i- Taking the Laplace transform in ti,. . . ,tn and rearranging the order of integration 
yields 

/■oo /"OO , 

C/(Ai,...,A„;mi,...,m„) = Af / .../ s^'-' . . . s^"-^Hi^){{s,, X,}^^i)dsi, . . . ,ds^. (18) 

Jo Jo 



Now multiply through by (t>{Xi + • • • + An) and apply Theorem 2.1 to obtain 



(f>{Xi H h A„);7(Ai, . . . , A„; mi, . . . , m„) (19) 

(^A + ... + _|_^^({5^,A.}r=i)dsi,...,ds„. 



6 



poo poo / poo 

-mY] ... ( / 

,_i Jo Jo \Jo 



s^'-K..s^--^-^H(^){{s,,X,}^^,)ds, ] dSi...ds,^ids,+i...dSn. (20) 



Let us now focus on the inner- most integral above. If nii = 1, we have 



dsi 



n 



d_ 


1 



H^"K{s^,K}7=l)ds^ 



i/("-i)(K,Afe}fe^,). (21) 



Above we have used the fact that exponentially as Si — > cx) , which follows from (|5|, and that 

H^K{s, = 0}) = i-i?(^)({sfc,Afe}fe^,) from [a). 



If > 1, we integrate by parts and get 



sT^-\..s::--'^Hi'-){{s,,X,yU)ds, = Y[s^^-' I sr'-'^HCn){{s,,X,}-^,)ds 



ds 



k^i 



,-1 9 



ds. 



= -("^.-1)114"'""' / sr-'HM{{s,,X,}-^,)ds, 

k^^ 



(22) 



With (21 1 and (22 1, (18) implies that (19 1 and (20 1 can now be written as 

0(Ai H h A„)J7(Ai, . . . , A„;toi, . . . ,m„) = ^miC/(Ai, . . . , A„;mi, . . . ,mi_i,TOi - l,mi+i, . . . , 

which finishes the proof. 



rrin), 



Remark. Theorem 3.1 is equivalent to Theorem 2.1 in [TS]. There, the result gives the joint moments of 
increments of the inverse subordinator: En"=i(-^(^i) ~ ^(•Si))'^S where < si < < • • • < s„ < t„, and 
ki, i — I . . .n are positive integers, and expresses them with an integral expression. The formulas given 
here, which are obtained using different methods, give directly Ei?(ti)™i . . .E{tn)™'" and have the added 
advantage of expressing higher order moments in terms of lower order ones. Also, out results give the Laplace 
transform of all n-time moments. 



Recursion: Theorem |3.1| gives expressions for the Laplace transform of all n-time moments of a general 
inverse subordinator. In theory, these n-dimensional Laplace transforms can be inverted, however doing so 
directly can be a formidable task. As an alternative, we take advantage of the recursive nature of equation 
(17). Let TV S Z be the order of the moment U{ti, . . . ,tm mi, . . . , m„), i.e. N — mi + mi -I- • • • -I- m„. 
Observe that equation ( 17l implies that moments of order N can be calculated using a linear combination of 
convolutions involving moments of order iV — 1 and the inverse Laplace transform of the function l/c/), which 
we address below. Thus, if one has C/(i) = V(t\ 1), then all moments can be obtained inductively using this 
method. 
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3.1 Getting the inverse Laplace transform 



To calculate moments of inverse subordinators, we must be able to invert Laplace transforms of the form f /4>, 
where the function / is known. Hence, we must first obtain the inverse Laplace transform of the function 
Following [7], section IILl, define the renewal measure to be the Borel measure whose distribution 
function is given by U{t). From this we see that for a.e. < a < 6, 

l(a,6](r)d[/(r) - U{b)-U{a) 

^ / P[E{b) > s] - P[E{a) > s]ds 
Jo 

= / P[D{s) <b]- P[D{s) < a\ds 
Jo 

' / V ( — 00,6] 

Jo 

l(a,6](^(s))rfs. 



E 



= E 



(23) 



By approximating with step functions, this can be extended to g{T)dU {t) — E Jq°° g{D{s))ds where g is 
a continuous function. Choosing g{T) — e^^^, we get that the Laplace transform of the renewal measure is 
given by 

Jo Jo 4>W 



Thus, for an arbitrary function /, the inverse Laplace transform of f /(j) is given by the convolution of / 
with the renewal measure, i.e. 



1 



/(A) 



f{t-T)dU{T). 



(24) 



This can be generalized to n dimensions. For ease of notation, write dU (t) — U'{T)dT, where U' is interpreted 
as a generalized function(since U might contain jumps). If < ti < t2 < ■ ■ ■ < tn, then 



(t)i\i + A2 + • • • + A„) 



/(Ai, A2, . . . , A„) 



'1 

"'0 

ti rt2 



f{tl - Ti, ■ ■ ■ ,tn - Tn)£ 



(/)(Ai + • • • + A„) 



{ti, . . . , tn) 

{Ti,...,Tn)dTi,...dTn 



JO 



f{tl - Tl, . . . , i„ - Tn)U'{Tl)S{T2 - Tl) . . . S{t„ - Tl)dTldT2 . . . dr^- 



This follows from standard rules of Laplace transforms (see for example [TH], page 169). This can be 
understood by writing C'^ — C'^^C'^^^_^ . . -J^^l then computing C'^^[\/(j){\i + • • • + A„)] by viewing 
A2 + • • • + A„ as a shift. Using the fact that ti < t2 • • • < in, the above reduces to 



c- 



1 

0(Ai + A2 + • • • + A„) 



/(Ai,A2 



1 A„) 



{ti,...,tn) — 



f{h-TM-r,...,U,-T)dU{T). (25) 
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Using /({A;}) = U{{Xi;mi}) from Theorem 3.1 equation (25 1 lets one compute any order moments by 
successive convolutions. As an application, we give a general expression for the covariance of an inverse 
subordinator in terms of U and dU. 



Corollary 3.2 Let D be a Levy subordinator with Levy exponent (j), and let E be the inverse subordinator 
of D. For s,t > 0, the covariance r, of E is given by 

/•sM 

Cov{E{s),E{t)) = {U{s-T) + U{t-T))dU{T)-U{s)U{t). (26) 
Jo 

Proof. We have 

GoY{E{s),E{t)) = ¥.E{s)E{t)~¥.E{s)¥.E{t) (27) 
= U{s,t- 1,1) ~U{s)U{t), (28) 

where we have used the notation ([T6|. From Theorem |3.1[ the Laplace transform of U{s,t; 1, 1) is 
C/(Ai,A2;l,l) = ^^^^-^^(c/(Ai,A2;l,0) + L/(Ai,A2;0,l)) 

0(Ai + Aa) y X2 ^1 / ' 

Now, ( |25[ ) gives 

/■sAt 

C/(s,t;l,l)= / iU{s-T) + U{t-T))dU{T). 







This with (28 1 establishes (26 1 



Remark. Let {D{s), s > 0} be a strictly increasing Levy subordinator with inverse {E{t),t > 0}. Then 
using Corollary |3.2[ one can show that E has stationary and uncorrelated increments if and only if D{s) = Cs 
for some constant C > 0. 



4 Examples 

Examples are discussed in |19j . We focus there on the following three important families of subordinators: 
• Poisson and Compound Poisson processes. Their Levy exponent is given by 

/•OO 

(/)(A) mA+ / {l-e-^'')i^{dx), 



where j/ is a finite measure on [0,oo). 

"Mixed" a-stable processes. Their Levy exponent is given by 

0(A) = C X^dpiP). 







Here p is some probability measure on (0,1). Notice the a-stable subordinator corresponds to the 
choice p{f3) = 5{a — j3). 
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• Generalized Inverse Gaussian Levy processes, which are Levy processes whose time 1 distribution is 
given by the Generalized inverse Gaussian distribution. This family of distributions was shown to be 
infinitely divisible in [3] and its Levy-Khintchine representation is derived in |10j . This family includes 
the gamma process, inverse Gaussian process and the reciprocal gamma process. 

In some cases, an analytic expression for the mean first-hitting time U(t) — EE(t) can be given for the 
above processes, however, this is usually not the case. In [19], we give numerical methods for inverting the 
Laplace transform U{\) = (A0(A))~^. We test these methods in cases where U can be computed explicitly. 
We then study the three examples above in detail, focusing in each case on the asymptotic behavior of U 
and on computing U and higher order moments numerically. 

A Appendix: Path Properties of Inverse Subordinators 

We describe here the path of the inverse subordinator. Figure [l] illustrates the relationship between a 
subordinator, {D{s), s > 0}, and its inverse, {E{t), t > 0}. Observe that both {D{s), s > 0} and {E{t), 
t>Q\ arc right continuous. 





















7. 







Figure 1: An example of a sample path by a subordinator D together with its inverse E. 

The following proposition provides additional details. We provide a proof for the convenience of the 
reader. 



Proposition A.l The sample paths of the inverse subordinator E are non- decreasing and are right contin- 
uous with left limits. The sample paths of E moreover, are continuous if and only if D is strictly increasing 



Proof. We start by proving that D increasing implies that E is non-decreasing and has cadlag sample 
paths (right continuous with left limits). Indeed, if ti < t2, then {s : D{s) > ^2} C {s : D{s) > ti}, meaning 
E{ti) < E{t2)- To prove the sample paths are cadlag, note that the left limits follow from the fact that E 
is non-decreasing. To see right continuity, observe that 

{s:D{s)>t}^ [j{s:D{s)>T}, 

^Sometimes, the inverse is defined as inf{s : D{s) > 0}. This process will have all of the same properties as the inverse 
defined here, writh the exception that it will be left continuous as opposed to right continuous. 
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since D is non-decreasing. Thus, taking inf's, 



E{t) = inf{s : D{s) > t} 

S 

= m{[j{s : D{s) > t} 

T>t 

= inf inf{s : D{s) > r} since inf{s : D{s) > r} G M for each t, 
— hni E{t), since inf{s : D{s) > rj increases in t. 

Thus, hm^^j+ E(t) — E(t), implying right continuity. 

To prove the second part of the proposition, we will first show that if E is continuous, then D is strictly 
increasing. Suppose the opposite is true, that is, D is not strictly increasing. If D{t) ~ c for ti < t < t2 with 
ti < t2, then by the right continuity of D, E{c—) < ti and E{c+) > t2, contradicting the continuity of E. 

We now show that D strictly increasing implies that E is continuous. Suppose that this is not that case, 
that E is discontinuous at some point c. From above, E is continuous from the right, thus discontinuity 
implies that E is not left-continuous at c, meaning E{c) — E{c—) > 0. Then, for any t with E{c—) < t < 
E{c) = inf{s : D{s) > c}, one has D{t) < c. Since E is non-decreasing, for c' < c one has 

t > E{c~) > E{c') = inf{.s : D{s) > c'}, 

thus D{t) > c'. Letting c' ^ c proves that D{t) = c, contradicting that D is strictly increasing. This finishes 
the proof. 



The next proposition displays the inverse relationship between the Levy subordinator and its first passage 
time. 

Proposition A. 2 Fix si, S2, ■ ■ ■ , Sn > 0. Then 

P[D{s,) <U, i = l,...,n] =P[E{U) > s,, z= for a.e. h,...,tn > 0. (29) 



Moreover, if D is strictly increasing, then {29) holds for all ti, . . . ,tn > 0. 

Proof. First, observe that we have the following set inclusions: 

{D{s,) <t„ ^= l,...,n} C {E{t,) > s,, i^l,...,n}(l {D{s,) < t,, i ^ l,...,n}. (30) 

To see this, suppose that D{si) < ti. Then by right continuity, D{s) < ti for s > Si sufliciently close to Si. 
Thus, E{ti) > Si since D is non-decreasing. 

For the second inclusion, assume E{ti) > s,;. If D{si) > ti, then we would have E{ti) — inf{s : D{s) > 



ti} < Si, a contradiction. This verifies (30 1 



Thus, we have for all si, S2, ■ ■ ■ , Sn > and ti, . . . , t„ > 0, we have 

P[D{s,)<t„ i^l,...,n]<P[E{t,)>s„ t^l,...,n]<P[D{s,)<t„ i^l,...,n]. (31) 
Now, for si, S2, ■ ■ ■ , Sn > fixed, 

n 

P[D{s^) <t„ 1^1,..., n] < P[D{s,) <t,, i^l,...,n] + J2 P[D{s.,) = t,]. (32) 
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For each i, let Ai = {t : P[D{si) — t] > 0}. Notice that each Ai is at most a countable subset of M+, and 
hence A = {{ti, . . . , t„) : ti ^ Ai for some i} is a set of Lebesgue measure 0. Thus, (32 1 implies that 

P[D{s,) <U, z = l,...,7i] < P[D{s,) < U, i = !,...,«], for (ti,...,t„) e A" 



This combined with (31) establishes (29 1 for all (ii, . . . G A"^. 

To prove the second statement of the proposition, assume D is strictly increasing. With this, the first 
inclusion in (30 1 is strengthened to equality: 

{D{si)<U, i^l,...,n} = {E{t,)> s,, i = l,...,n}. (33) 

To see this, we proceed by contra-positive. If D{si) > ti, then since D is strictly increasing, D{s) > ti for all 
s > Si, hence E{ti) — inf{s : D{s) > ti} < Si. This, combined with (30) proves (33), thus P[D{si) < ti, i — 
1, . . . ,n] = P[E{ti) > Si, i = 1, . . . ,n] for all ii, . . . , t„ > 0. 
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